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The usual mathematical formulation of availability 
assumes an exponential distribution for failure and repair 
times. While this is sometimes correct for reliability, it 
is not likely to be for maintainability. This study was 
conducted no verify that the lognormal ana gamma distribu- 
tions are suitable descriptors for corrective maintenance 
repair times, and to estimate the difference caused in 
assuming an exponential distribution for availability and 
maintainability calculations when in fact the distribution 
is lognormal. Forty- six sets of data of electronic and 
mechanical systems and equipments were analyzed using the 
methods of probability plotting and statistical testing for 
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I • INTRODUCTION 

a. OVERVIEW 

The effectiveness of a system depends not only on its 
ability to meet its specified performance requirements, but 
also on its ability to perform when needed, for the duration 
of its assigned missions, and for its operational lifetime. 
The technical disciplines concerned with these time-related 
system characteristics are reliability, maintainability and 
logistic support. 

Commonly used methods for prediction and demonstration 
of corrective maintenance repair times assume the validity 
of the lognormal distribution. 

In order to help in focusing attention on these matters, 
a statistical analysis on data sets of demonstrated and 
field repair times has been conducted as part of a study on 
the application of the lognormal and gamma distributions to 
corrective maintenance data. The results of the analysis 
are given in this thesis. A previous study of electronic 
systems and equipments concluded that the lognormal dis-ri- 
bution was a suitable descriptor for repair times (Ref. 1]. 
B. PURFOSE AND APPROACH 

1 . Ob iecti ves 

The objectives of this study have been; 
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(a) To test tha exponential, lognormal, and gamma 
distributions as descriptors of repair rime data on mechani- 
cal and electronic systems and equipments, 

(b) to verify that the lognormal or gamma distri - 
bution is a suitable descriptor for repair times, 

(c) to estimate the percentage error caused in 
assuming exponentiali ty for availability and maintainability 
calculations when in fact the distribution is lognormal. 

2. Systems and D at a Analyzed 

Forthy-six sets of repair time data for electronic 
and mechanical systems / equipments were analyzed (Table 1) . 
Data for the electronic systems (Sets 8-20) came from formal 
maintanability demonstration test reports furnished to us by 
the Rome Air Development Center, U.S. Air Force. Mechanical 
equipments included field repair time pumps used in French 
nuclear electric power stations (Sets 1-7), and additional 
26 sets of field repair data from British Aircraft, provided 
by the United Kingdom Ministry of Defence, representing 
elapsed times (Sets 21-33) , and man hours (Sets 34-46). 'Jo 
detailed reports are available for the mechanical items. 

3 • Analysis Ap proa ch 

Three approaches were used for testing data. 
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(a) Rough check using a technique called Boxplot. 
This gives a rough idea that the data comes 
from symmetric (normal) or skewed distributions 
lognormal, exponential or gamma). 

(b) Plotting the data using appropriate probability 
paper. 

(c) Statistical tests of significance for assumed 
distributions. 

The following procedure was used for analyzing the 

dat a, 

i. Sketch boxplot and histogram of data and 
determine if data appears to come from a symmetric or skewed 
distribution. Obtain an idea about center and dispersion of 
dat a. 

ii. Plot the data on appropriate probability 
paper for the exponential and lognormal distributions. In 
some cases, plots were mads for the gamma distribution. 

iii. A Chi-Square Go cdness -0 f-Fit test was per - 
formed for lognormal and exponential assumptions, [Ref. 2] 

iv. A W test, due to Shapiro and Milk, [Ref. 3] 
was used to test the lognormal assumption for sample sizes 
less than 51, (due to availability of the tables). 
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v. A Kolmogoro v-Smirnov (K-S) test was used for 
lognormal and gamma assumptions [Ref. 4]. 
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II -STATISTICAL CONSIDERATIONS 
A. ROUGH CHECK METHODS 

1. Box pl ot (Refs. 5 and 6] 

Purpose: Many batches of data pile up in the middle. 
To analyse the behavior of a batch, we need a picture of 
where the middle (median) lies and how -he tails relate to 
it. The middle is generally better defined than the nails 
so we want to see more at the tails. Some values called 
outliers are so low or so high that they seem to stand apart 
from the rest of the batch. 

One technique for displaying such batches of data 
which is often mere convenient than a histogram is the box- 
plot. But the histogram is still good for showing the shape 
of the distribution. Boxplot will not show multiple modes. 
There may be several "modes”, and we want to see them. 

A boxplot is obtained by first calculating the lower 
and upper quartiles and the median of the batch (sample) of 
numbers and then plotting th a numbers on a horizontal line. 

The lower quartile is the value that divides the 
batch into two parts, with 1/4 of the numbers below this 
value and 3/4 above it. The upper qaartile is the value 
with 3/4 of the observations below and 1/4 observations 
above it. 
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The next step is to draw a narrow rectangular box 
with ends corresponding to the lower and upper guartiles, 
and to display the median point by a plus sign. The length 
of this box is called the interquartile range or H-Spread 
(Figure 1) . 

Tc the above are added lines marking 1/2 H-Spread on 
each side of the box. Dana values outside the lines (out- 
liers) are marked with stars and rectangular circles. 

*** 1 + I * * 0000 

interquartile range 
\.5 interquartile * range 

i interquartile range ' 

-Figure 1 - Boxplot- 

First, the data are ordered such that, 

X(1) < X (2) < X (n) where, 

X(1) : smallest (minimum) value 

X(n) : largest (maximum) value 
median®!! , center of batch 
if n is an odd number, M=X((a + 1)/2) 

if n is an even number , M= (X ( (N/2) + 1) +X ( (n/2) - 1) ) /2 
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guartiles=H , center of half of batch 
d (H) = (| d (M) | ♦ 1) /2 where | A| means integer part of A 
inter guartile distance = upper guartile -lower guartile 

Here are good boxplot examples for the exponential 
and lognormal model (data generated randomly). 



BOXPLOT EXPONENTIAL 

-I + I 0 

-Figure 2 - Boxplot Expor er.tial- 
BOXPLOT NORMAL (log is taken then boxplot) 

X + x 

-Figure 3 - Boxplot Normal- 



2. Hi stogram [Ref* 5 and 6] 

Boxplot is good for comparison but a histogram is 
still good for showing shape. Sometimes histogram is more 
helpful than the boxplot, especially for detection of more 
bumps in the data. Histogram prints data in classes. Exam- 
ples are given below in Figure 4. 

HISTOGRAM EXPONENTIAL 



MIDDLE OF 
INTERVAL 

0. 

1. 

2 . 

3. 

4 . 

5. 

6 . 

7. 

8 . 



NUMBER OF 
OBSERVATIONS 

g ********* 

1 0 ** ** * **** * 

2 ** 

0 
0 
0 
0 

1 * 
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HISTOGRAM NORMAL (log is taken then histogram) 



MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERVATIONS 


1.0 


1 


* 


1.5 


3 


*** 


2.0 


6 


***** * 


2.5 


6 


***** * 


3.0 


5 


***** 


3.5 


2 


** 


4.0 


1 


* 


4.5 


2 


** 



-Figure 4 - Histograms- 



B. STATISTICAL DISTRIBUTIONS 



1 . The E xp onenti al Dis tribut ion 



The probability density function of the exponential 



distribution is (Ref. 2], 



f(x,^)=Ae , x > 0 and A> 0 

Its cumulative distribution function is, 

? (X ,\) = \ ^dt = 1-e 

which can be easily evaluated. 



The exponential distribution is used frequently as a 



time to failure model for a system when a constant failure 



rate is assumed [Ref. 4]. But for repair times, it can be 



shown that it is not an appropriate model. Because repair 



time includes diagnostic, correction, and verification tasks 



a repair time density muse have a value of zero at time t =0 . 



It should then increase to its maximum value rapidly and 



then gradually decrease towards zero as time increases. 
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The maximum likelihood estimator of the parameter of 
:he exponential distribution from given data is, 



a 1 

A = where 



X X(i) 
i. — 1 



x=- 



2. Lognormal Distribution 

This distribution has many different shapes for non- 
negative variates. It is skewed to the right, the degree of 
skewness increasing with increasing values of standard devi- 
ation. The mean and standard deviation are scale and shape 
parameters respectively. 

A random variable is said to have a lognormal dis - 
tribution if the logarithm of the variable is normally dis- 



tributed, 

1 

f (x) ■== 

xG/rff 

^ 0 

with mean and variance 

2 

fl = E (lnX) an d (J = 




Varian 



e(lnX) 



If data X (1 ) , X ( 2) , 



X (n) and if 



Y ( 1) =ln X ( 1) , Y (2) =1 nX (2) , Y(n)=lnX(n) than, 

maximum likelihood estimators of the lognormal distribution 
are, 

A a2 1 n 2 

JLl= y C= — £ (y(i)-y) 

n i = 1 
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3 



th 



The p 


percentile is calculated 


by [ Refs . 7 


Median 


Z ( . 5 ) =0 X 

0.5 


e ^ 


th 


Z (. 9) =1 .28 2 X 


£* + 1. 282 (7 


90 per. 


e 


th 


0.9 


+ 645^ 


95 per. 


Z (.9 5) = 1.645 X 

0.95 



* Ga mma Distribution 

A random variable x has a gamma distribution, and it 
is referred to as a gamma random variable, if and only if 
its probability density is given by, [Ref. 9] 

1 a -1 -x/b 

e for x > 0 

elsewhere 




where a > 0 and b > 0 

a is shape parameter 
b is scale parameter 
Pis the gamma function 

The mean and variance of the gamma distribution are 
given by. 



2 

mean=a.b and variance=a.b 



Special cases of the gamma distribution play an im- 
portant roles in statistics. For instance, for a=1 and b=Q, 
we get the exponential distribution with X=1/Q. 
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met ers 



For the gamma distribution the scale and shape para- 
can be estimated from data in many ways; here are two 
a. Method Of Moment (M.O.M) estimators 



2 

(sample mean) 

a = — - 

sample variance 

where mean=a.b 



b = 



sample variance 
sample mean 
2 



and variance=a. b 



b. Maximum Likelihood Method (M.L.S.) estimators 
This method is more complex as compared to the 
M.O.M.; it requires numerical calculation. There is an APL 
function named ESTGAM that gives directly M.L.E. (App.C). 

It is used because it provides a better estimator of the pa- 
rameters a and b in large samples, provided the data truely 
come from a gamma model. 

C. PB03ABILITY PLOTTING 

Plotting data points on probability paper is quite easy 
and does not require many calculations. It is a visual exa- 
mination rather than statistical calculation (Figure 5) . 

What can we learn from a probability plot? 

1. Visual indication of assumed distribution. 

2. Examine departures from assumed distribution. 
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5a. Tails lonaer than normal 
5b. Tails shorter than normal 

-Figure 5- 

1 • Lognormal Distribution 

Consider the order statistics of a sample 

x(1) < x(2) < < x (N) Then, x (i) 

is an estimate of the (i/N) th percentile of X 
If X is N (0,1) then we expect 
-1 

x(i )2zJ?r (i/N) for i = 1 , 2 , 3 N 

-1 

Then a plot of X (i) vs jo (i/N) should give a 
straight line. The only problem is 

-1 -1 

Jcf (i/N) ^infinity, so we use if (i/(N+s)), 3 = 1. 
PROCEDURES 

a. Use normal probability paper 
Plot (i/N+s) vs X (i) 

b. Use r.ormplot in APL workspace [Ref. 21] 
Use HINITAB 

N SC ORE S of Cl, put in C2 

PLOT C2 vs Cl 

NSCORES are (i/N + s) 
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Consider log of data Y < 1 ) < Y (2) < 

< Y (N) where Y(i)=LnX(i) and then apply above 

procedure to get lognormal plot. 

2. Sx^cn en tial Distribution 

Suppose X is exponential A with sample size n, so 
Pr ( X < x) = 1-exp (-A X) where X > 0. Then (i/r.+s) is 

an estimate of Pr (X < X ( i) ) =1-exp (-Ax (i) ) 

So if the data really come from an exponential dis- 
tribution (i/N+s) vs. 1- e xp (-Ax (i) ) , should yield a 
straight line. 

PROCEDURES 

a. In (i/N+s) vs. X (i) on rectangular paper, 

b. EXPON PLOT in APL workspace, [Ref. 21] 

c. MINITAB, 

d. Use Gamma probability paper for ALPH A = 1 
BET A=Q where Q is exponential parameter so that, 

f (x) = (1/Q) *exp (-X/3) for X > 0 

3. Gamma Distribution 

There is an APL program function named 3AM that uses 
estimated gamma parameters and gives test statistic D. The 
program also gives a matrix for probability plotting pur- 
poses. (see Appendix C) 
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D. STATISTICAL TESTS FOR DISTRIBUTIONAL ASSUMPTIONS 

We need to know whether we are actually dealing with a 
sample from an exponential population or whether the data 
values of random variables have a lognormal or some ocher 
distribution. A statistical rest of significance of a dis- 
tributional assumption provides an objective technique for 
tentatively assessing whether or not an assumed model pro- 
vides an adequate description of the observed data. 

There are three basic steps involved in statistical tes 
methods . 

1. A test statistic is calculated from the data, 

2. The probability of obtaining the calculated test 
staris-ic is determined, 

3. Assessment is made of the adequacy of the assumed 
distr ibution. 

a. If the probability of obtaining the calculated 
test statistic is "LOW'*, one can say that the assumed dis - 
tribution dees not provide an adequate representation. 

b. If the probability of obtaining the calculated 
test statistic '•HIGH" , then the data provides no evidence 
that the assumed distribution is not adequate, judged by th 

Z. 93 U • 
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The definition of "LOW" or "HIGH" depends on the 
user’s preferences, and the consequences of rejecting the 
distribution. For example, the rejection statistic is typi- 
cally be 0.01, 0.05, or 0.1. For the purpose of this study 
0.05 was selected as the reject criterion. 

It should be pointed cut that a statistical test, 
although it allows one to reject an assumption as inadequate 
does not prove that the assumed distribution is correct. 

1 • Th“ Chi.- Squared Goqdness-of - Fit Test 

This test is one of the oldest and most commonly 
used for evaluating distributional assumptions. Basically, 
the given data are grouped into cells and compared to the 
expected number of observation in the cells based on the 
assumed distribution. Then if the test statistic, calcula- 
ted from this comparison exceeds a critical chi-square value 
the assumed distribution is rejected. 

The only problem is dividing the data into cells. 

If the number of observations is small it is usually sugges- 
ted to use a number of cells as large as possible, subject 
to the restriction that it should not exceed N/5 ( N is sam- 

ple size). 
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The computations involved in the Chi-Squared Good - 
ness-Of-Fit test were made by a comparer program, in which 
are given as inputs the assumed cumulative distribution 
function, the number of observations, the number of equipro 
bable cells, and the number of parameters estimated from th 
sample. The outputs are the Chi-Square statistic and its 
probability of exceeding a chi-square value for a given nun 
ber of degrees of freedom. 

The procedure used is described as follows : * Ref . 2] 
a. The cell boundaries are determined from the 
assumed cumulative distribution as the values such that the 
estimated probability of the observation value falling 
within a given class is 1/K. For each class, 

Pr (X<x (i) ) = i/k, where 

x=the random observation to be assigned to the 

i-th cell 

x (i) =the i-th cell boundary to be solved from 
the above formula 

k = the number of cells, in this study k = r./5 

The lower bound of the first c=ll and the upper 
bound of the last cell are the smallest and the largest 
values of the observations. 
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b. 

cell E (i) 

c. 

is counted 

d. 



The expected number of observations for each 
is equal to n/k for each cell. 

The number of observed values in each cell M(i) 
based on the results of the above equation. 

The test statistic is 



x^-i i 



i= 1 



. 2 T 
M(r) I 



n 



e. The computed value )( is 
level of significance, or the probabili 
value with degrees of freedom equal to 
number of estimated parameters from the 
calculated chi-square 



used to compute the 
ty of a 3hi-Square 
K-M-1 where H is "he 
sample to exceed 



2 2 

<*=Pr( (V) > X ) with V degrees of freedom 

1 -pc 

If ck is less than or equal to .05 the assumed dis- 
tribution can be rejected as inadequate. 

2* The Kolmoqoro v- Smi rnov Test (K-S Test) 



K-S test examines the cumulative distr ibut ion func- 
tion rather than probability density function . It simply 
compares the sample cumulative distribution function S(.) 
to the hypothesized distribution function F(.). As a mea- 
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sure of comparison the test uses, D = Max | F (. ) -S ( . ) i ; this 
is the main difference between the Chi-Squared 3.0.?. test 
and the K-S. test. [Ref. 4] 

On the other hand here we are dealing with Type I 
error, that is the hypothesis is true and reject it. The 
probability of commiting a Type I error is called the sigr.i 
ficance level for the test alpha. Here alpha is taken to b 
0.05. This means no more than a 5 % chance of rejecting th 
true hyphothesis. 

X has CDF F (x) = Pr ( X < x) , 

Take a sample X(1) < X(2) < < X(N), 

We want to estimate F(x)= Pr(X < x) = S (x) 

n 



S (x) 


= 0 


for X < 


X (1 ) 






n 

If 


1/N 


" X > 


X ( 1 ) . AND. X 


< 


X(2) 


It 

«« 


2/N 

• 


»• X > 


X (2) . AND. X 


< 


X (3) 


It 

If 

ft 


K/N 

• 


" X > 


X (K) . AND. X 


< 


x (K + 1) 


ft 

ft 


• 

1 


X > 


X (N) 






Find 


D= MAX| 


1 S (X) -F (X) 1 








Using 


Table 


9 of reference 4 for gi 


.vsn sign 



lev el , 

If D exceeds the tabulated value that comes from th 
table, reject the hypothetical distribution. 
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If D does not exceed the tabulated value chat comes 
from table, do not re jeon the hypothetical distribution. 

There are two kinds of tables for the K-5 test. One 
of these is used if the parameters are estimated from data, 
and the other one is used if the parameters are given. Here 
parameters are estimated from data. 

3. The W test for lognormal issu motion 

The ^ test is shown in Ref. 3 to be an effective 
procedure for evaluating the assumption of normality against 
nonnormal alternatives, even if only a relatively small num- 
ber of observations are available. Hahn and Shapiro suggest 
that the w test may also be used to evaluate the assumptions 
of the lognormal distribution. This follows because of the 
property that if the logarithms of the observations follow a 
normal distribution, then the original values of the obser - 
vations are loqnormally distributed. 

The following procedure was aDpliea using tables 
from Reference 2. 



a. The observations are ordered. 



X (1) < X (2) < . . . < X (i) 



< X (n — 1 ) < X (n) 



b. The following parameters are computed. 



i 



n 




.£ lnx(i) 
1=1 



n 
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ii. If n is even, k=n/2 

n is odd, k= (n- 1) /2 



B= a (k-i+ 1) ,[ (Inx (n-i+1) ) - 

where the values of a(k-i+1) fo 

are given in table IX [Ref. 2] for n = 3 to r. = 50 . 

iii. The test statistic W is 

2 2 
W= B / S 

c. The approximate probability of obrai 
calculated value of W can be obtained from the f 



Z=A+^ln 



!:§! 

1-w| 



using the values of A ,TJ and 6 given in Table XI 
for the appropriate sample size and then using s 
normal distribution to determine the probability 
ing a value less than or equal to Z, which is th 
cance level of the test. 

q = Pr(Z < Z) 

g 

If q is less than or equal to .05, the s 
vel of significance in this study, the lognormal 
as an inadequate assumption. 



(Inx (i) ) ] ] 
r i= 1 , 2 , . . k 



ning the 
or mu la 

of Ref. 2 
tar.dardized 
of obtain- 
e signifi- 

elected le- 
is rejected 
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III. DAT A A: MUSIS RESULTS 

A. RESULT OF TESTS FOR DISTRIBUTIONAL ASSUMPTIONS 
1 . S ummary of re su its 

Tables 2,3 and 4 summarize the results of the sta - 
tistical test analysis. In all cases the 5 % level was cho- 
sen as the reject criterion. The 46 sets of daza used in 
-his study can be divided into three parts as indicated in 
Section on IB2. 

The following is an example of the computer program 
summary table which includes zhe results for the Chi-Square 
test and calculated parameters from the sample data. 

SET NO 18 - AN/ GXS-2 (V) 



SAMPLE SIZE 


N = 26 


NO. OF CELLS 


K = 5 


(a) 


SAMPLE MEAN 


= 19.70 


STANDARD DEV 


= 22.62 


(b) 


EXPONENTIAL 


LOGNORMAL 


DIFFERENC 


u 


PARAM1 


0.05 


2.56 




(C) 


PARAM2 




0.78 






MTTR 


19.70 


19. 13 


2.9 6 * 


(3) 


50-TH PERCNT 


13.65 


12.93 


5.61 % 




90-TH PERCNT 


45.3 5 


40. 22 


12.77 % 




95-TH PERCNT 


59.00 


55.46 


6.39 % 




CHI-SQR STAT 


5. 9 2 


0 . 54 




(e> 


DEG OF FREED 


3 


2 




(f) 


SIGNIF LEVEL 


0. 115E+00 


0. 764S+00 




(g) 
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Table 2 



TOTAL DATA ANALYSIS RESULTS 

Set Sample EXP. LOG NORMAL GAMMA 



NO 


Size 


TEST 






TEST 










TEST 






Pr. 


Ch. 




Bo. 


Pr. 


Ch. 


W K 


-S 




M. 


M. 






PI. 


Sq. 


R e j . 


PI. 


PI. 


Sq. 


m 

X • 


T. 


Rej. 


OM 


LS 


1 


43 


3 


R 


+ 


G 


G 


A 


A 


a 




R 


R 


2 


43 


B 


R 


+ 


G 


J 


A 


R 


A 




R 


R 


3 


30 


B 


S 


+ 


G 


F 


A 


A 


R 




R 


R 


4 


43 


B 


R 


+ 


B 


B 


R 


V 


R 


+ 


R 


R 


5 


33 


3 


R 


+ 


E 


? 


E 


R 


A 


% 


R 


R 


6 


21 


B 


A 




B 


3 


R 


R 


R 


+ 


R 


A 


7 


59 


3 


R 


+ 


P 


F 


R 


— 


R 


+ 


R 


R 


8 


20 


G 


R 


+ 


P 


p 

4. 


R 


R 


R 


+ 


A 


A 


9 


49 


D 


R 


+ 


G 


a 


A 


A 


R 




A 


A 


10 


50 


P 


R 


+ 


F 


F 


R 


A 


a 


* 


R 


R 


11 


50 


P 


R 


+ 


G 


n 


A 


A 


A 




A 


A 


12 


43 


3 


R 


+ 


F 


F 


R 


A 


R 


* 


A 


A 


13 


30 


B 


R 


+ 


P 


p 


R 


R 


a 


+ 


A 


A 


14 


26 


P 


R 




G 


a 


A 


A 


A 




A 


A 


15a 


33 


G 


R 


+ 


P 


r* 

a 


A 


A 


R 




R 


R 


16 


30 


B 


R 


+ 


B 


3 


R 


R 


R 


+ 


R 


R 


17 


50 


P 


R 


+ 


G 


j 


A 


A 


A 




R 


A 


18 


26 


F 


A 




F 


F 


A 


A 


A 




A 


A 


19 


57 


B 


R 




B 


B 


A 




R 


+ 


A 


A 


20 


57 


B 


R 


+ 


G 


F 


R 




R 


+ 


E 


A 


21 


157 


P 


R 


+ 


G 


F 


R 




R 




E 


R 


22 


28 


F 


A 




G 


J 


A 




A 




R 


A 


23 


36 


B 


R 


+ 


P 


B 


R 


R 


A 


* 


P 


a 


24 


58 


B 


R 


+ 


B 


a 


R 




R 


+ 


R 


R 


25 


20 


F 


R 


+ 


G 


B 


R 


“a 


A 




A 


A 


26 


92 


F 


R 


+ 


F 


9 


R 




R 


+ 


R 


R 


27 


21 


B 


R 


+ 


B 


B 


R 




a 


* 


A 


A 


28 


33 


B 


R 


+ 


G 


jr 


A 


A 


A 




c 


A 


29 


47 


P 


R 


+ 


G 


p 


A 


A 


A 




R 


A 


30 


35 


3 


R 


+ 


F 


a 


A 


A 


R 




A 


A 


31 


81 


P 


R 


+ 


B 


3 


R 




a 


+ 


A 


A 


32 


22 


B 


R 


+ 


G 


o 


A 


HC 


R 




A 


A 


33 


45 


B 


R 


+ 


B 


B 


R 


R 


a 


* 


T? 


a 


34 


171 


3 


R 


+ 


G 


F 


R 




R 


+ 


R 


R 


35 


31 


F 


A 




G 


G 


A 




A 




A 


A 


36 


37 


3 


R 


+ 


P 


J 


A 


A 


A 




R 


R 


37 


60 


■p 


R 


+ 


3 


B 


R 




a 


+ 




a 


38 


20 


G 


R 


+ 


G 


F 


A 




A 




A 


A 


39 


95 


B 


R 


+ 


F 


G 


A 




A 




R 


R 


40 


22 


T5 


A 




G 


■j 


A 


“X 


A 




A 


A 


41 


35 


F 


R 


+ 


G 


G 


A 


A 


A 




■o 


A 


42 


51 


F 


A 




G 


ST 


A 




A 




r\ 

A 


A 


43 


38 


B 


R 


+ 


F 


F 


A 


"TT 


R 


* 


A 


A 


44 


87 


F 


R 


+ 


p 


J 


A 




a 




A 


A 


45 


22 


F 


R 


+ 


G 


G 


A 


~X 


a 




A 


A 


46 


47 


B 


a 


+ 


B 


3 


R 


R 


R 


+ 


R 


R 



2 

Criteria for rejection Pr (X ) or Pr (W) less than .05 

or D value for K-S test bigger than tabulated value. 
* To be determined from probability olot or SoxDlot 
where B=2ad, P=?oor, 7=Fair, G=Good 

K = Do not reject, R=Reject 
Pr . P 1 . = Probability plot, Bo.Pl.= Soxolot 
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T a ble 3 

LOGNORMAL AND EXPONENTIAL RESULTS 





rSampIe 


EX'?. 


rourr 


ThFl 






p 


Size 


2 


2 




1 


K-S 


Test. 




P(X ) 


P(X ) 


D 


(W) 

“07535 “ 


Ts. St. 


Tab. 7a. 


i 


53 


— ururrg -- 


“* -07659" 


— 


7TIT<d 


7T35 


2 


43 


0. 00166 


0. 358 




0.025 


.07 9 


. 135 


3 


30 


0. 48E-4 


0. 0588 




0.0745 


.205 


.161 


4 


43 


0.59E-7 


0.0005 




0.0 


. 238 


. 135 


5 


33 


0.018 


0. 0128 




0. 0328 


.148 


. 154 


6 


21 


0.245 


0.0468 




0.0002 


.27 


. 193 


7 


59 


0.238E-6 


0.596E- 


-7 




.172 


.115 


8 


20 


0.010 1 


0.0143 




0.017 


.235 


. 190 


9 


49 


0.1 19E-6 


0. 459 




0.31 


.217 


.126 


10 


50 


0.5 96 E-7 


0. 0445 




0.099 


.212 


. 125 


11 


50 


0.000 11 


0.979 




0.632 


.0 86 


.125 


12 


43 


0.596E-6 


0. 0421 




0. 19 


.143 


. 135 


13 


30 


0. 28E-5 


0. 01 86 




0.01 


.3 


.161 


14 


26 


0. 129E-4 


0.292 




0. 238 


. 151 


. 174 


15 


50 


0.208 E-4 


0. 00408 


0.001 


.245 


.125 


15a 


33 


0. 2 16 E-4 


0. 484 




0. 06 


.211 


.154 


16 


30 


0. 8 11 E-4 


0. 032 1 




0.0 


.199 


.161 


17 


50 


0.0323 


0.494 




0. 583 


. 1 25 


. 125 


18 


26 


0.115 


0.76 4 




0.36 


. 084 


.174 


19 


57 


0 . 596 E-7 


0. 360 






. 23 


. 1 17 


20 


57 


0.0 


0. 116E- 


-4 




.44 1 


. 117 


21 


157 


0.0 


0.0 






.073 


.07 


22 


28 


0. 381 


0.52 6 




0.624 


.145 


.167 


23 


36 


0. 808E-4 


0. 173S- 


-3 


0. 0074 


. 132 


. 147 


24 


58 


0. 132E-3 


0.02 5 5 






.143 


.116 


25 


20 


0. 0334 


0.91 22- 


-2 


0. 35 


.16 


.190 


25 


92 


0.6 56 E-6 


0.380E- 


-4 




.102 


.092 


27 


21 


0. 1 74E-2 


0.231E- 


-2 


0. 175 


. 30 1 


. 193 


28 


33 


0.018 


0. 554 




0.345 


. 143 


. 1 54 


29 


47 


0.69E-2 


0.71 1 




0. 524 


.106 


. 130 


30 


35 


0.254 E-4 


0.592 




0.6 17 


.263 


.149 


31 


81 


0. 186 E-3 


0.86 IE* 


-4 




.121 


.098 


32 


22 


0.468E-2 


0. 06 3 




0.59 


.219 


.189 


33 


45 


0.4 27 E-3 


0.315E- 


-3 


0.013 


. 147 


. 132 


34 


17 1 


0.5 96 E-7 


0. 596E- 


-7 




.069 


.068 


35 


31 


0.456 


0.724 




0.97 


.117 


. 159 


36 


37 


0.0258 


0.57 1 




0.11 


.083 


.145 


37 


60 


0.542E-5 


0. 430E- 


-2 




. 1 23 


.114 


38 


20 


0.0334 


0.094 




0. 896 


.099 


. 110 


39 


95 


0. 298E-5 


0. 22 






.091 


.091 


40 


22 


0.071 6 


0.201 




0.424 


.156 


.189 


41 


35 


0.0476 


0. 126 




0.48 


. 125 


. 149 


42 


51 


0.193 


0. 103 






.092 


. 1 24 


43 


38 


0.470E-2 


0.24 2 




0.04 


.265 


. 144 


44 


87 


0. 182E-2 


0. 193 






. 106 


. 095 


45 


22 


0.463E-2 


0.340 




0.57 


.224 


. 189 


46 


47 


0.7 57 S-5 


0. 192E- 


-3 


0. 004 


.179 


.13 
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19 
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21 

22 

23 

24 
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I able 4 

GAMMA DISTRIBUTION TEST RESULTS 
(Kolmogorov Smirnov G.O.F.Tast) 



D value of 


A/R 


D value of 


A/3 


Tabulated 


K. 0. M . Sst. 


_ __ 


M . L . E . Est. . 




Value 


. 36 2 


3 


. 154 


a 


.135 


.286 


R 


. 182 


R 


.135 


. 22 


R 


. 185 


3 


.161 


.45 8 


R 


. 344 


R 


.135 


.37 


R 


. 21 


R 


. 1 54 


. 23 


R 


.164 


A 


.193 


. 29 9 


3 


. 138 


R 


.115 


. 18 


A 


. 154 


A 


.190 


.09 3 


A 


. 1 1 


A 


.126 


. 143 


R 


. 147 


R 


.125 


. 08 1 


A 


. 08 1 


A 


.125 


. 124 


A 


. 107 


A 


. 135 


.148 


A 


. 152 


A 


.161 


. 168 


A 


. 153 


A 


. 174 


. 158 


R 


. 162 


R 


.125 


. 17 


R 


.114 


R 


. 154 


. 23 1 


R 


. 162 


R 


.161 


. 136 


R 


. 097 


A 


.125 


. 23 2 


R 


. 143 


A 


.174 


. 073 


A 


. 074 


A 


. 1 17 


. 145 


R 


. 109 


A 


. 117 


. 192 


R 


. 163 


R 


.07 


. 12 


A 


.112 


A 


.157 


. 285 


R 


. 20 1 


3 


.147 


.26 


E 


. 164 


3 


.116 


. 172 


A 


. 155 


A 


.190 


.22 2 


R 


. 14 


R 


.092 


. 186 


A 


. 187 


A 


. 193 


. 172 


R 


. 09 


A 


.154 


. 187 


R 


.113 


A 


. 1 30 


. 1 


A 


. 093 


A 


.149 


.036 


A 


. 085 


A 


.098 


. 105 


A 


.114 


A 


.189 


. 185 


R 


. 142 


R 


. 132 


.284 


R 


.13 


D 

_l 


.068 


. 153 


A 


. 098 


A 


. 159 


. 263 


R 


. 158 


R 


.145 


. 286 


R 


. 172 


a 


.114 


. 124 


A 


. 107 


A 


.19 


.336 


R 


. 161 


3 


.091 


. 08 


A 


.085 


A 


. 1 84 


. 154 


R 


. 099 


A 


. 149 


. 163 


R 


. 099 


A 


.124 


. 13 1 


A 


. 126 


A 


. 144 


. 08 8 


A 


. 075 


A 


.095 


. 12 


A 


. 08 9 


A 


. 189 


. 23 9 


R 


. 175 


R 


.13 
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Notes : 



a. The number of equiprobable cells, X, was chosen 
as N/5 where N the sample size. 

b. The sample mean and standard deviation are cal- 
culated based on the maximum likelihood estimates. 

c. PA3AH1 for the exponential distribution is the 

reciprocal of the sample mean, for the lognormal distribu- 

A. AO 

tion PAR AH 1 and PAR AH 2 are pi and U , the parameters of the 
lognormal distribution. 

d. The HTTR and the 50th, 90th and 95th percentiles 
are calculated from the sample and are based on the rela- 
tionships between the calculated parameters and their dis- 
tribution functions. The percentage difference is between 
the exponential and lognormal HTTR and percentiles. 

e. The Chi-Square statistic. 

f. The number of degrees of freedom (K-H-1) is K-2 
for the exponential distribution and K-3 for the lognormal. 
This is because one parameter is estimated from the sample 
in the exponential case and two in the lognormal case. 

g. The level of significance is the probability of 
a Chi-Square variate with the specified degrees of freedom 
exceeding the calculated Chi-square statistic. 
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For the exponential model, only the Chi-Squared test 



was used. In 6 out of 46 oases we would accept the exponen- 
tial model. If we compare the exponential vs. lognormal mo- 
del using the Chi-Square test, test prefers the lognormal 
model tc the exponential model. In only one case would the 
exponential model be accepted and the lognormal rejec-ed. 

ex ponen t ial 

A R 

A (5 |19| 24 

lognormal I ~ T -- I ~ ~ ~~ I 

R | 1 | 2 1 | 22 

6 4 0 

Therefore the lognormal model is a better descriptor 
with respect to the exponential model. 

In most cases the lognormal assumption works reason- 
ably well. In those cases in which the results of the chi- 
squared test, W test and K-S test do not agree, a plot on 
lognormal probability paper was used to determine the appro- 
priateness of this assumption (* in Table 2). Comparison 
between test methods were made in three ways for the lognor- 
mal model and tabled below. 

a. Chi-Sq. vs. K-S test (46 cases) 

K-S 

A R 



A 


I 15 | 


1 10 | 


25 


Sq. 

R 


t — 1 
1 3 ! 


1 

1 18 | 


21 



18 2 3 
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b. Chi-Sq. vs. W test (34 crises) 

W 



test methods sometimes 
and the chi-square 
chi-square agrees with 
the W test in 82 % of the cases. H and K-S test agree 62 %. 
Therefore, there is no significant difference among the 
three test methods. 

A previous study has been dene using 24 data sets of 
electronic systems and equipments, using only the W and chi- 
squared tests for the lognormal model [3ef. 1]. The cable 
below gives a comparison between the exponential and lognor- 
mal for the previous study. For the exponential model do 
not reject (A) or reject (3) decision was made cased on the 
chi-squared test. For the lognormal model accept by either 
the chi-square or W test was used as a criterion. 



Chi- Sq. 



19 | 2| 

1 



21 

13 



JLi— ?-i 

23 1 1 

c. W vs. K-S test (34 cases) 



K-S 

A R 



A 

R 



23 

11 



13 | 10 | 

16 1 8 

For the lognormal methods the 
give different results. The K-S test 
test agree in 74 % of the cases. The 
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ex ponen tial 
A R 



lognor mal 



I 7 

!;t 

8 



14 | 

;r| 

1 6 



21 

3 



The rabies below show same kind of results 
current study. Only difference between the two st 
there are different kind of systems end equipments 
mechanic or electronic) in the current study. 



logncr mal 



A 

R 



ex ponential 


A 


R 


1 5 


23 


1 1 


1 7 


6 


4 0 



29 

18 



For all cases (46 cases) 



ex ponent ial 
R 



lognormal 



A 

R 



I 1 
1 _ 0 . 
1 



7 I 

:r! 

1 2 



8 

5 



For electronic equipments only 



These results indicate that the locrnormal 
better than the exponential model for repair time 
For the gamma modal, the ilLS estimator is 
werful than the MOM estimator, with the K-S test u 
testing the data. The decision was made to accept 
suits of the MLE when the two estimators gave diff 



for the 
udies is 
(e.g. 



(13 cases) 

model is 
data, 
more do - 
sed for 
the re - 

6T 
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results, since in no case did the MOM suggest acceptance 
when the MLE said to reject. In the eight cases where there 
was a difference, therefore, the MLE was used as the crite - 
rion. This gives the following results for the gamma model. 



MLE 

Therefore using the MLE, we would accept the gamma 
model in 27 out of 46 cases. In 16 of these 27 cases, the 
lognormal model is also accepted and in 5 additional cases 
all three distribution assumptions would be accepted. Again 
only one case (case 6) would we accept the gamma and expo- 
nential model but reject the lognormal. 

We can also compare lognormal and gamma model using 
the K — S test. 



MOM 



1 

I 

1 

r*U | 


a 


19 ! 

I 


l 

00 I 

1 

1 


0 1 


1 19 



19 27 



lognor mal 



gamma 
A a 



A 


13 1 


1 5 

t 


R 


13 | 


15 




26 


2 0 



13 

29 



The above table shows than the K-S rest prefers the 
gamma model. 27 out of 46 cases accepted the gamma model, 
but only 13 out of 46 cases accepted the lognormal model. 
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Tables 6-9, (Appendix B) summarize the results for 
all cases and for each of the three different types of data. 
Table 10 (Appendix B) gives total results in terns cf sample 
sizes of 20-29, 30-39, 40-50 and over 50. For exponential 
assumption, increasing sample size increases rejection per - 
centages. For gamma assumption, increasing sample size also 
increases rejection. In lognormal, increasing the sample 
size decreases rejection up to sample sizes of forthy, then 
rejection increases rapidly. 

Because of the large number of data sets, only some 
interesting cases of the data analysis and probability plots 
for those sets are discussed in the next section. 

2. Di scuss ion of sped fic cases 

The following discussions illustrate some cases 
which are of special interest. The data analysis results 
and plots are given in Appendix A. 

a. Set No. 18 - AN /GXS- 2 (V) 

This case is an illustration where all three 
models are accepted (not rejected) for all tests. Data were 
collected from 26 maintenance tasks. In this case, as in 
most other of the formal maintainability demonstration (ca- 
ses 8-20) , the tests were performed in accordance with the 
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U.3. military standard for maintainability demonstration 
(MIL- STD-47 1) [Ref. 23] in which sample selection was ge - 
nerated from knowledge or estimates of component failure 
rates and tests were conducted using service technicians of 
the appropriate skill levels and given prior training on the 
equipments . 

This case fits the lognormal, exponential and 
also gamma for both estimators. All three tests give high 
significance for lcgncrmali tv. 

b. Set No. 6 - Reactor Coolant Pump 

This case is the only one the exponential model 
accepted and the lognormal model rejected. However the 
plots are not good. Further investigation is necessary in 
this case but since only the data was available, it was not 
possible for us to do so. 

c. Sets No. 28 and 41 - Main Rotor Head and Blades 
These two cases indicate a very good example of 

logncrmalit y . All three tests signal "do not reject" fcr 
the lcgrcrmal assumption. The plots are also very good. 

The gamma family does not well represent the data sets. 

d. Sets No. 1 and 2 - French Fessenheim Pumps 
( Down / Repair ) times 
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These -wo cases are discassed as unsuccessful 
for gamma and exponential but good for lognormal since at 
least 2 out of 3 tests give acceptance for the lognormal 
assumption. For repair times the w test rejects the assump 
tion while the Chi-squared and K-S tests do not, but down 
times are rejected by only the K-S test. There are some 
outliers around 5%, but the rest of the data show a good 
lognormal fit in the probability plot. 

The existence of outliers could be a reason for 
rejection by one of the 3 tests. For downtime data, after 
trimming the two high outliers the lognormal assumption was 
also accepted by the K-S test. For repair time, after trim 
ming one outlier the data shows high significance (0.15) fo 
the W test. Thus, in some cases it would be valuable fo 
examine in more detail the reason for the high or low oufli 
ers . In this case, only the data was made available so tha 
this was not possible. 

e. Set No. 4 - Heater Drain Pump 
Set No. 5 - Main Feedwater Pump 
Set No. 4 illustrates the case where all three 
tests for the lognormal model are rejected. Here outliers 
are important factors. The effect of outliers can be seen 



41 



easily from the boxplot, histogram and probability plot. No 
5 is almost the same as no 4, the only difference between 
the two being that the K-S test barely accepts the lognor - 
mality in set no 5. For these cases, trimming the outliers 
did not change anything in so far as significance levels or 
acceptance or rejection are concerned. 

f. Set No. 9 - AN/ARC-162 (Intermediate) 

Set No. 10 - AN/ARC-163 (Intermediate) 

These two cases are discussed together, because 
they give opposite results. Derailed reports exist in both 
cases. Test technicians were used for running the tests and 
removina the faults. [Refs. 12 and 13] 

The exponential assumption is rejected in bc^h 
cases. In set No 9 the lognormal assumption is rejected only 
by the K-S test while the other tests give high signifi- 
cance. The probability plot and boxplot also signal that 
'•do not reject” is appropriate for this case. The gamma 
model fits this set for both estimators. 

In Set No. 10 only the S test accepts, but the 
Chi-Squared test barely rejects (0.045) and the K-S test 
rejects the lognormal model. However the boxplot and proba- 
bility plot are fairly good for this set. An examination of 
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the report indicates that some bias exists by the use of 
'•average” times for some repair time segments for much of 
the data instead of normal data. Therefore, a further exa - 
mination of this data vould be in order, 
g. Set No. 15 - AN /UYK- 1 4 ( C) 

Set No. 15a - -Revised AN/UYK-14 (C) 

Set No. 15 illustrates what happens when the da- 
ta points are not selected at random and thus introduce bi - 
as. In this case, 33 tests were selected according to the 
sample selection procedure in MIL-STD-471A [Ref. 23]. How- 
ever in order to have 50 data points according to the test 
plan 17 additional data points using ''average” repair times 
for several components from a previous test were used. This 
therefore introduces a bias in the data, resulting in the 
initial rejection of the lognormal model {Table 3). After 
taking out the biased data points the analysis shows that 
the lognormal model is a good descriptor. The histogram 
does not look very good, but the probability does for the 
lognormal model. 

This example then shows that if data does not 
fit an assumed distribution and if the reason is bias, then 
after finding a good reason for throwing out some bias do - 
ints, the data frequently will agree with a lognormal model. 
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h. Sets No 21 and 34 - Fuselage Structure 

These two cases are example of large sample si - 
zes (157 and 171). For both data sets all distributional 
assumptions were rejected by all test methods. Bur for the 
lognormal case probability plots, boxplots and histograms 
are fairly good. These cases should be further invesziga - 
ted . 



B. E SHOES IN CALCULATED PARAMETERS 

1 . Error in MTTR when Assu ming. an E xpone ntial 
Distribution Inste a d of Lognormal 

The percent error in the mean-time-to-repair (MTTR) 
is calculated as follows. 



M 

log 



t he 



E 



M - M 
log exp 



M 

lo g 



lognormal mean = s 



x 100 
fi + O.SG* 



where , 



M - the exponential mean = x 
exp 

The results of the percentage error in the mean 
which are summarized in Table 5. All cases have an error 
less than 18 X. Therefore since the inherent availability 
formula normally used is. 



MTBF 1 

A = = 

I MTBF + MTTR MTTR 

1 + 

MTBF 
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T able 5 

PERCENTAGE ERROR IN MTTR WHEN ASSUMING AN 
EXPONENTIAL INSTEAD OF LOGNORMAL DISTRIBUTION 



Set 

No 

-1- 


Sample 

Mean 

(EX?) 

-2- 


Sa mple 
Mean 
(LOG) 
-3- 


Error (%) 
13-21 
3 


1 


119. 51 


112. 4 4 


6.29 


2 


24. 0 


22. 1 7 


8. 25 


3 


121. 63 


140. 34 


13.63 


5 


57. 21 


48.6 9 


17.49 


9 


7. 23 


7. 22 


0. 16 


10 


9. 35 


9.34 


0.08 


1 1 


43. 05 


4 3.59 


1 . 23 


12 


15. 75 


15.74 


0.02 


14 


15.73 


15. 54 


1 .25 


15 


15. 85 


15.59 


1.65 


17 


17. 42 


17.76 


1.93 


18 


19. 7 


19.13 


2.96 


22 


1. 87 


2.0 7 


9.73 


23 


1. 55 


1.4 8 


4.9 


25 


1. 76 


1.7 8 


1.36 


27 


0. 57 


0.3 1 


1.36 


28 


0. 41 


0. 28 


1.29 


29 


2. 45 


2.44 


0.69 


30 


1. 7 


1. 7 2 


1 .45 


32 


1.64 


1.66 


0.85 


35 


4.03 


4.57 


11.84 


36 


2. 55 


2.3 9 


6.96 


38 


2. 26 


2.3 3 


2.96 


39 


4. 22 


4. 1 


2.73 


40 


2. 5 


2. 57 


2.57 


4 1 


3.5 


3.5 1 


0.27 


42 


3. 63 


3. 66 


0. 63 


43 


2. 4 


2.4 4 


1.78 


44 


1. 72 


1.37 


7.88 


45 


2. 58 


2.60 


0.58 



NOTE: Only cases included when locrnormal model was 

accepted ( at least one out of three test methods ) 
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and since for a reasonable availability MTTR « MTBF (e.g. 
MTTR is measured in hours, MTBF is 100*s of hours), then it 
is shown that the exponential assumption for STIR introduces 
negligible error in availability. 

C. CONCLUSIONS 

From the data analysis in this study as well as the pre- 
vious study [Ref. 1], it is concluded that the lognormal 
model is a good descriptor for repair times. 28 out of 46 
sens show that, within a 0. 05 level of significance for at 
least one of the 3 test methods , this model can not be re- 
jected. Similarly, the data analysis shows that the assump- 
tion of the exponential distribution should be rejected in 
40 sets. The gamma assumption is not rejected in 27 cases. 

The percent difference in the MTIR, when assuming expo- 
nential distribution when the true distribution is lognormal 
has been found to be negligible for calculated system avai- 
lab il ity. 

When the results of the statistical tests indicate oppo- 
site conclusions, probability plots, boxplots and histograms 
are helpful in determining an acceptable model and to explo- 
re anomalies in the data. Probability plots are also useful 
in estimating percentiles, as well as density parameters. 
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SET NO 1 


APPENDIX A 

- FRENCH FESSENHEI M POMPS 


(DOWN TIME) 


SAMPLE SIZE 


N = 43 


NO. OF CELLS 


K = 7 


SAMPLE MEAN 


= 1 19. 51 


STANDARD DEV 


= 225.79 


EXPONENTIAL 


LOGNORMAL 


DIFFERENCE 


PARAM1 


0.0 1 


3.89 




PARAM2 

MTTR 


119.51 


1 . 66 
112.44 


6.29 % 


50 -TH PERCNT 


82.84 


49. 10 


68.70 % 


90-TH PERCNT 


275.1 9 


255.75 


7.6 0 % 


95-TH PERCNT 


358.0 2 


408.08 


12.27 % 


CHI-SQR STAT 


13.16 


2.42 




DEG OF FREED 


5 


4 




SIGNIF LEVEL 


0. 219 E- 01 


0.659E+00 




30X PLOT SET 1 

1+ I * 


O 0 




0 



30XPL0T LN SET 1 
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HISTOGRAa SET1 



MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERV ATIONS 


10.00 


6 


***** * 


15.00 


6 


***** * 


20. 00 


0 




25.00 


1 


* 


30.00 


3 


*** 


35.00 


2 


** 


40. 00 


1 


* 


45.00 


0 




50.00 


0 




55.00 


4 


**** 


60.00 


4 


** ** 


65.00 


0 




70.00 


0 




75. 00 


0 




80.00 


2 


** 


TOGRAM LNSET1 




MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERV ATIONS 


1.5 


1 


* 


2.0 


0 




2.5 


1 1 


***** **** * * 


3.0 


2 


** 


3.5 


6 


***** * 


4.0 


8 


** *** ** * 


4.5 


3 


** * 


5.0 


6 


***** * 


5.5 


2 


** 


6.0 


2 


** 


6.5 


1 


* 


7.0 


1 


* 
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SET NO 2 


- FRENCH FESS 


ENHEIM PUMPS 


(REPAIR TIME) 


SAMPLE SIZE 


N = 43 


NO. OF CELLS 


K = 7 


SAMPLE MEAN 


= 24.00 


STANDARD DEV 


= 31.95 


EX EON ENT I AL 


LOGNORMAL 


DIFFERENCE 


PARAM1 

PARAM2 


0.0 4 


2.73 

0.73 




MTTR 


24. 00 


22. 17 


8.25 % 


50-TH PSRCNT 


16.6 4 


15.40 


7.99 % 


90-TH PERCNT 


55. 26 


46.00 


20.14 % 


95-TH PERCNT 


71.90 


62.70 


14.67 % 


CHI-SQR STAT 


19.35 


4.37 




DEG OF FREED 


5 


4 




SIGNIF LEVEL 


0. 166 E-02 


0. 358E+00 




PLOT SET2 


-I + I 


* O 


O 


C 0 


PLOT LN SET 2 


± 


+ I 
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HISTOGRAM SET 2 



MIDDLE OF 
INTERVAL 
0 . 
20 . 
40. 
60. 
80. 
100 . 
120 . 
140. 
160. 



NUM3ER OF 
OBSERV ATI ON S 



12 

24 

3 

1 

1 

0 

0 

1 

1 



************ 

**************** ******** 
** * 

* 

* 



* 

* 



HISTOGRAM LNSST2 



MIDDLE OF 
INTERVAL 

1.5 
2.0 

2.5 

3.0 

3.5 

4.0 

4.5 

5.0 



NUMBEP OF 
OBSERV ATI ONS 

U **** 

9 ******** 

1 1 *********** 

1 i ***** ****** 

5 ***** 

1 * 

1 * 

2 ** 
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HePAitt TJmGS 



2*i 



iOo' 



.5^ KJ o L 

Xoa^omnol Plot 



Ac 



$ 

*••• 



•• 



• • 
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• • 



» • 
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SET NO 4 



HEATER DRAIN POMP 



SAMPLE SIZE 


N = 43 


NO. OF CELLS 


K = 7 


SAMPLE MEAN 


= 67.60 


STANDARD DEV 


= 154.16 




EXPONENTIAL 


L03N0RMAL 


DIFFERENCE 


PARAM1 




0.0 1 


2.90 




PARAM2 






1.81 




MTTR 




67.5 0 


44.72 


51.18 % 


50-TH PERCNT 




46. 86 


18.09 


> 100 % 


90-TH PERCNT 




155.6 7 


101.51 


53.3 6 % 


95-TH PERCNT 




202.53 


165.42 


22.43 % 


CHI-SQR STAT 




56. 4 7 


20.00 




DEG OF FREED 




5 


4 




SIGNIF LEVEL 




0. 596 E-07 


0.500E-03 




BOXPLOT SET4 










+ IOO 


0 




0 0 0 


0 


BOXPLOT IN SET 4 










r 




+ I 


* * 


0000 
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HISTOGRAM SET4 



MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERVATIONS 


5.00 


9 


***** **** 


10. 00 


12 


***** ******* 


15.00 


7 


***** ** 


20.00 


6 


***** * 


25. 00 


2 


** 


HISTOGRAM LNSET 


2 




MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERVATIONS 


1.5 


5 


***** 


2.0 


9 


***** **** 


2.5 


9 


********* 


3.0 


12 


** ** * ***** ** 


3.5 


1 


* 


4.0 


1 


* 


4.5 


0 




5.0 


2 


** 


5.5 


0 




6.0 


2 


** 


6.5 


2 


** 
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SET NO 5 



MAIN FEEDWATER PUMP 



SAMPLE SIZE 


N = 33 


NO. OF CELLS 


K = 6 


SAMPLE MEAN 


= 57.21 


STANDARD DEV 


= 105.37 


EXPONENTIAL 


LOS NORMAL 


DIFFERS 


PARAM1 


0.0 2 


3.20 




PARAM2 




1 . 38 




MTTR 


57. 2 1 


48.69 


17.49 % 


50-TH P5RCNT 


39. 6 6 


24 . 47 


62.03 % 


90-TH PERCNT 


131.74 


110.10 


19.65 % 


95-TH PERCNT 


171. 39 


168.54 


1.6 9 % 


CHI-SQR STAT 


11.91 


10.82 




DEG OF FREED 


4 


3 




SIGNIF LEVEL 


0. 180 E- 01 


0. 128E-01 




BOX PLOT SST5 








1+ I — 


O O 


O 




BOX PLOT LNSET5 








T 


+ I 
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HISTOGRAM SET5 



MIDDLE 0? 


NUMBER 


OF 


INTERVAL 


OBSERVA 


lTIONS 


5.00 


2 


** 


10.00 


10 


***** **** * 


15.00 


3 


*** 


20.00 


3 


*** 


25.00 


0 




30.00 


5 


** *** 


35.00 


1 


* 


40. 00 


1 


4c 


45.00 


2 


4C 4c 


50.00 


1 


4c 


55.00 


0 




60.00 


1 


* 


TOGRAM LNSET 


5 




MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OESERV ATI ONS 


1.0 


1 


4c 


1.5 


1 


4c 


2.0 


2 


4c* 


2.5 


10 


***** **** * 


3.0 


4 


**** 


3.5 


7 


***** ** 


4.0 


4 


** ** 


4.5 


0 




5.0 


0 




5.5 


2 


** 


6.0 


2 


** 
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IMPAIR TINS (4 tor 4 ) 
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SET NO 6 



REACTOR COOLANT PUMP 



SAMPLE SIZE 


N = 21 


NO. OF CELLS 


II 

& 


SAMPLE MEAN 


= 53.57 


STANDARD DEV 


= 39.22 


EXPONENT IAL 


LOGNORMAL 


DIFFERENCE 


PARAM1 


0.0 2 


3.64 




PARAM2 




1 . 49 




MTT 3 


58.57 


80.08 


26.86 % 


50-TH PERCNT 


40.6 0 


38.05 


6.69 % 


90 -TH PERCNT 


134. 97 


181.78 


25.81 % 


95-TH PERCNT 


175.46 


283.05 


38.01 % 


CHI-SQR STAT 


2. 8 1 


3.95 




DEG OF FREED 


2 


1 




SIGNIF LEVEL 


0. 245 E+00 


0.468E-01 





BOX PLOT SET6 



-I + 



30XPL0T LNSET6 



* * * I + I 
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HISTOGRAM SET6 



MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERV ATI ONS 


0 . 


4 


**** 


20. 


0 




40. 


3 


*** 


6 0. 


6 


***** * 


30. 


3 


** * 


100. 


4 


**** 


12C. 


0 




140. 


0 




160. 


1 


* 


TOGRAM LNSET6 




MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERV ATI ONS 


0.5 


1 


* 


1.0 


0 




1.5 


2 


** 


2.0 


1 


* 


2.5 


0 




3.0 


0 




3.5 


3 


*** 


4.0 


7 


***** ** 


4.5 


6 


***** * 


5.0 


1 


* 
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CwMMtAr<V« Pl^tfNT AH 



SET NO 9 - AN/ARC-162 (V) (INTERMEDIATE LEVEL) 



SAMPLE SIZE 


N = 49 


NO. OF CELLS 


K = 10 


SAMPLE MEAN 


= 7. 23 


STANDARD DEV 


= 3.44 


EXPONENTIAL 


LOGNORMAL 


DIFFERENCE 


PARAM1 


0. 1 4 


1 . 38 




PARA M2 




0. 19 




MTTR 


7.23 


7. 22 


0.16 % 


50-TH PESCNT 


5. 0 1 


6.57 


23.64 % 


90-TH PERCNT 


16.65 


11.49 


44.93 % 


95-TH PERCNT 


21.6 7 


13.46 


60.99 % 


CHI-SQR STAT 


47.94 


6.71 




DEG OF FREED 


3 


7 




SIGNIF LEVEL 


0. 119 E-06 


3. 4593+00 





BOX PLOT SET9 



I + I ** * 



BOX PLOT LN SET 9' 



I + I 
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HISTOGRAM SET9 



MIDDLE 0? 
INTERVAL 

2 . 

4. 

6 . 

8 . 

10 . 

12 . 

14 . 

16. 

18 . 



NUMBER OF 
OBSERV ATI ON S 
0 

12 ************ 

18 ****************** 

5 ****** 

5 ***** * 

3 *** 

1 * 

2 ** 

1 * 



HISTOGRAM LN5ET9 



MIDDLE OF 
INTERVAL 
1.0 
1.2 

1.4 
1 .6 
1.8 
2.0 
2.2 

2.4 
2.6 
2.8 



NUMBER OF 
OBSERV ATION3 
2 ** 

3 *** 

5 ***** 

5 ***** 

14 ************** 

7 ******* 

4 ** ** 

4 ** ** 

2 ** 

3 *** 
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CWMDkATlVC P«^t*MTACk 



SET NO 10 - AN/ ARC- 16 3 (7) (INTERMEDIATE LEVEL) 



SAMPLE SIZE 


N = 50 


NO. OF CELLS 


II 

o 


SAMPLE MEAN 


= 9. 35 


STANDARD DEV 


= 4.27 


EX PON ENT IAL 


LOSNORMAL 


DIFFERENCE 


P ARAM1 


0. 1 1 


2.14 




PARA M2 




0. 18 




MTTR 


9. 35 


9.34 


0.08 % 


50 -TH PERCNT 


6. 4 8 


8.53 


24.03 Z 


90-TH PERCNT 


21.52 


14.73 


46.09 % 


95 -TH PERCNT 


23.00 


17.20 


62.92 % 


CHI-SQR STAT 


53. 2 0 


14.40 




DEG OF FREED 


8 


7 




SIGNIF LEVEL 


0. 596 E- 07 


0. 445E-01 




PLOT SET 10 








T 


+ I 




** 



BOX PLOT LN SET 10 



I + I 
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HISTOGRAM SET 10 



MIDDLE OF 
INTERVAL 

4. 

6 . 

8 . 

10 . 

12 . 

14. 

16. 

18. 

20 . 

HISTOGRAM LNSET 

MIDDLE OF 
INTERVAL 
1.2 

1.4 
1.6 
1.8 
2.0 
2.2 

2.4 
2.6 
2.8 
3.0 



OF 

ATI ON S 
***** 

***** **** * 

** ** * *** * ** ** * * 
********* 

** 

* 

** * 

*** 

** 



OF 

ATI ONS 
* 

***** *** 

** ** * ** 

***** ******** ** 
**** 

* 

***** 

*** 



NUMBER 

OBSERV 

5 

10 

15 

9 

2 

1 

3 

3 

2 

10 

NUM3ER 

OESSRV 

1 

3 

3 

9 

7 

15 

4 

1 

5 

3 
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Rf»PAlf4> Tit*** (pr\lr%} 
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5 ST NO 15 - AN/ UYK- 1 4 (C) 



SAMPLE SIZE 


N = 50 


NO. OF CELLS 


o 

II 


SAMPLE MEAN 


= 15.85 


STANDARD DEV 


= 11.06 


EX PON ENT IAL 


LOGNORMAL 


DIFFERENCE 


PA RAMI 


0.06 


2.58 




PARAM2 




0. 33 




MTT 3 


15.85 


15.59 


1.6 5 % 


50-TH PERCNT 


10.99 


13.21 


16.85 % 


90-TH PERCNT 


36.50 


27.63 


32.09 % 


95-TH PERCNT 


47.48 


34.05 


39.45 % 


CHI-SQR STAT 


3 5. 5 0 


20.80 




DEG OF FREED 


8 


7 




SIGNIF LEVEL 


0. 208 E-04 


0. 408E-02 





BOXPLOT SET 1 5 



1 + I * * 



BOXPLOT INSET 15 



I + I 
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HISTOGRAM SET 15 



MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERVATIONS 


5. 


8 


***** *** 


10. 


1 8 


************ 


15. 


1 1 


***** ****** 


20. 


3 


*** 


25. 


3 


*** 


30. 


2 


** 


35. 


0 




40. 


1 


* 


45. 


4 


**** 


TOGRAK LNSET 


15 




MIDDLE 0? 


NUMBER 


OF 


INTERVAL 


OBSERVATIONS 


1.8 


3 


*** 


2.0 


5 


****** 


2.2 


12 


***** ******* 


2.4 


5 


***** 


2.6 


4 


**** 


2.8 


7 


***** ** 


3.0 


3 


*** 


3.2 


3 


*** 


3.4 


2 


** 


3.6 


1 


* 


3.8 


4 


**** 



74 







75 



SET NO 1 5 A - REVISED AN/3 YK- 1 4 (C) 



SAMPLE SIZE 


N = 33 


NO. OF CELLS 


K = 6 


SAMPLE MEAN 


= 14.37 


STANDARD DEV 


= 9.66 


EXPONENT IAL 


LOGNORMAL 


DIFFERENCE 


PA RAMI 


0.07 


2.50 




PARAM2 




0.31 




MTTR 


14. 37 


14.19 


1.29 % 


50-TH PE RC NT 


9.96 


12.15 


17.98 % 


90-TH PE SC NT 


33. 1 0 


24.83 


33.30 % 


95-TH PE3CNT 


43.06 


30.40 


41.65 % 


CHI-SQR STAT 


26.82 


2.45 




DEG OF FREED 


4 


3 




SIGNIF LEVEL 


0. 216 E-04 


0. 484S+00 





SOX PLOT SET 15 A 



-I + I * * * 0 

BOX PLOT LN SET 15 A 



1 + I 
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HISTOGRAM S ST 1 5 A 



MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERV ATI ON 5 


5. 


8 


***** * 


10. 


1 1 


***** * 


15. 


6 


** *** * 


20. 


2 


** 


25. 


2 




30. 


2 


** 


35. 


0 




40. 


1 


* 


45. 


1 


* 


TOGRAM LNSET 


15A 




MIDDLE OF 


NUM3ER 


OP 


INTERVAL 


OBSERV ATI ONS 


1.8 


3 


*** 


2.0 


6 


***** * 


2.2 


6 


** ** * * 


2.4 


4 


** ** 


2.6 


4 


** ** 


2.8 


2 


** 


3.0 


2 


** 


3.2 


2 


** 


3.4 


2 


** 


3.6 


1 


* 


3.8 


1 


* 
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Set Mt> 

_Lo^no«-*Tt«l Tlo-b 



t 



i 




s 

5 




X 9 to gm JO MO f» t# 7a 4a 5# x J§ 
CUMUcATivC ^*dcfHr>0< 
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SET NO 18 



AN/GXS-2 (V) 



SAMPLE SIZE 


N = 26 


NO. OF CELLS 


II 

cn 


SAMPLE MEAN 


= 19.70 


STANDARD DEV 


= 22.62 


EXPONENTIAL 


LOGNORMAL 


differs: 


P ARAM 1 


0.0 5 


2.56 




PARA M2 




0.78 




MTTR 


19.70 


19. 13 


2.96 % 


50-TH PERCNT 


1 3. 6 5 


12.93 


5.6 1 % 


90-TH PERCNT 


45. 3 5 


40.22 


12.77 % 


95-TH PERCNT 


59.00 


55.46 


6.3 9 


CHI-SQR ST AT 


5.92 


0.54 




DEG OF FREED 


3 


~> 




SIGNIF LEVEL 


0. 115 E+00 


0 . 7 64 E + 00 




BOXPLOT SET 18 








--I + I 




O 


O 


BOXPLOT LNSET 18 









I + I 
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HISTOGRAM S2T18 



MIDDLE OF 
INTERVAL 

0. 



NUMBER OF 
OBSERVATIONS 
2 ** 



10. 


14 


***** ********* 


20. 


5 


***** 


30. 


2 


** 


40. 


0 




50. 


0 




60. 


1 


* 


70. 


0 




80. 


0 




90. 


2 


** 



HISTOGRAM LNSET1S 



MIDDLE OF 
INTERVAL 
1.0 

1.5 
2.0 

2.5 

3.0 

3.5 

4.0 

4.5 



NUM3ER OF 
OBSERV ATIONS 
1 * 

3 *** 

g ****** 

6 ****** 

5 ***** 

2 ** 

1 * 

2 ** 
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ftcrAift tlMfiS (tn ifl.) 



Set Wo IB 

- } O^OOe-mal ?lo^ 



/* 



/ L-- _ , 

i 5 »• i* jo u » i« r» $9 i* 99 

CUHUUT|V£ PE« CENTAL 



31 



Set Vjo 46 
Gohnoaq "Plot 



n 

$503. 



•V 



5 

4 



Qwiu^ rawij. -aivrfsy 



32 



CuMUUAT/VtS pCfLcevvTACt. 



SET NO 21 



FUSELAGE STRUCTURE 



SAMPLE SIZE 


N = 1 57 


NO. OF CELLS 


K = 24 


SAMPLE MEAN 


= 1 . 54 


STANDARD DEV 


= 1.83 


EXPONENT IAL 


LOGNORMAL 


DIFFERENCE 


PARAM1 


0.6 5 


-0.06 




PARAM2 




0.97 




MTTR 


1.54 


1.53 


1.01 % 


50-TH P5RCNT 


1.07 


0.94 


13.52 % 


90-TH PERCNT 


3.5 6 


3.33 


6.93 % 


95-TH PERCNT 


4.6 3 


4.75 


2.64 % 


CHI-SQ3 STAT 


169.06 


155.0 




DEG OF FREED 


22 


21 




SIGNIF LEVEL 


0 . 0 


0.0 




PLOT SET21 

— I + I 




OOO 0 0 


0 00 



BOX PLOT LNSET21 



I + I 
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HISTOGRAM SET 21 

EACH * REPRESENTS 2 OBSERVATIONS 



MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERV ATI OHS 


0. 


33 


***** **** ****** ** 


1. 


73 


***** ** ************** *************** 


2. 


27 


** ** * ******** * 


3. 


7 


** ** 


4. 


4 


** 


5. 


2 


* 


6. 


5 


*** 


7. 


3 


** 


8. 


1 


* 


9. 


0 




10. 


2 


* 


TOGRAM LNSET 


21 




MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERV ATI ON S 


-2.5 


3 


*** 


-2.0 


2 


** 


- 1.5 


12 


** ** * ** ** * ** 


-1.0 


16 


** *** ********* * * 


-0.5 


34 


***** ** ** ************ ***** ******** 


0.0 


36 


***** ********** ************* ******** 


0.5 


27 


*************************** 


1.0 


10 


***** **** * 


1.5 


7 


***** ** 


2.0 


8 


***** *** 


2.5 


2 


** 
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SET NO 28 - 


MAIN ROTOR 


HEAD AND BLADES 


(ELAPSE TIH: 


SAMPLE SIZE 


N = 33 


NO. OF CELLS 


K = 6 


SAMPLE MEAN 


= 2. 07 


STANDARD DEV 


= 1.87 


EXPONENTIAL 


LOGNORMAL 


DIFFERENCE 


PAHA Ml 


0. 4 8 


0.46 




PAHAM2 




0.52 




MTTR 


2.07 


2.05 


0.89 % 


50-TH PE EC NT 


1.4 4 


1.58 


9.09 f 


90-TH PEHCNT 


4. 77 


4.00 


19.31 % 


95-TH PEHCNT 


6. 2 1 


5.20 


19.33 1, 


CHI-SQR STAT 


11.91 


2.09 




DEG 0? FREED 


4 


3 




SIGNI? LEVEL 


0. 180 E-01 


0. 554E+00 





BOX PLOT SET28 



— I + I 0 



BOX PLOT LNSET28 



I + I * 
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HISTOGRAM SET28 



MIDDLE OF 
INTERVAL 

0 . 

1 . 

2 . 

3. 

4. 

5. 

6 . 

7. 

8 . 

9. 



HISTOGRAM LNSET 

MIDDLE OF 
INTERVAL 
-0.3 
-0.4 
0.0 
0.4 
0.8 
1.2 
1.6 
2.0 
2.4 



OF 

ATI ON S 

***** *********** 
****** 

** *** ** 

** 

* 



* 



OF 

ATION5 

***** *** 
***** * 
***** * 

** ** * * 

* 

* 



NDM3ER 

OBSERV 

0 

16 

6 

7 

2 

1 

0 

0 

0 

0 

0 

1 

28 

NUMBER 

03SERV 

3 

2 

8 

6 

6 

6 

1 

0 

1 
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"iifcJPA.H. TlMJC 






Ho 2B 



• •• 



M 0 



b 



I _ ■ 

1 5 I# 29 S 9 M» u *# T 9 i O $9 $r J* 

COMuwAmvt. «m 
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CUMl/i-ATlVf *>6RcSNT“AGJi 



SET NO 34 



FUSELAGE STRUCTURE 



SAMPLE SIZE 


N = 1 71 


NO. OF CELLS 


K = 34 


SAMPLE MEAN 


= 2. 83 


STANDARD DEV 


= 4.60 


EXPONENTIAL 


LOGNORMAL 


DIFFERENCE 


PARAM1 


0.35 


0.36 




PARAM2 




1 . 26 




MTTR 


2. 83 


2.68 


5.26 % 


50-TH PERCNT 


1.96 


1.43 


36.77 % 


90-TH PERCNT 


6. 5 1 


6.03 


7.9 5 % 


95-TH PERCNT 


8.46 


9.05 


6.51 % 


CHI-SQR STAT 


126.6 5 


92.35 




DEG OF FREED 


32 


31 




SIGNIF LEVEL 


0. 596 E-07 


0.596E-07 





BOX PLOT SET34 



1+ I * **00 0 OOO 0 0 0 

BOX PLOT LNSET34 



* 1 + x ** 
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HISTOGRAM SET 34 



5 OBSERVATIONS ARE BELOW THE FIRST CLASS 



MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERVATIONS 


0. 500 


38 


***** **** ***************************** 


1. 000 


33 


*************************** ****** 


1. 500 


25 


***** ** ** **************** 


2. 000 


15 


***** ********** 


2. 500 


9 


***** ** ** 


3. 000 


1 3 


************* 


3. 500 


2 


** 


4. 000 


3 


*** 


4. 500 


3 


*** 


5. 000 


3 


** * 


5. 500 


0 




6. 000 


5 


***** 


6- 500 


1 


* 


7. 000 


1 


* 


HISTOGRAM LNSET34 




MIDDLE OF 


NUMBER 


OF 


INTERVAL 


OBSERVATIONS 


-2. 500 


2 


** 


-2.250 


0 




-2. 000 


0 




-1.750 


3 


*** 


- 1. 500 


5 


***** 


-1. 2 50 


0 




- 1. 000 


1 1 


***** ****** 


-0.750 


15 


***** ********** 


-0. 500 


6 


***** * 


-0. 250 


13 


***** ******** 


0. 0 


19 


** *** **** ****** ** * * 


0. 250 


15 


***** ********** 


0. 500 


19 


******************* 


0.750 


12 


***** ******* 


1. OCO 


18 


****************** 


1. 250 


3 


** * 


1. 500 


8 


***** *** 


1.750 


6 


***** * 


2.000 


2 


** 


2. 250 


7 


***** ** 


2. 500 


1 


* 


2.750 


3 


*** 


3.000 


0 




3. 250 


2 


** 


3.5 00 


1 


* 
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SET NO 41 - 


MAIN ROTOR 


HEAD AND BLADES 


(MAN HOURS) 


SAMPLE SIZE 


N = 35 


NO. OF CELLS 


K = 7 


SAMPLE MEAN 


= 3. 50 


STANDARD DEV = 


3.42 


EX PONENIIAL 


LOGNORMAL 


DIFFERENCE 


PARAM1 


.0.29 


0.90 




PARAM2 




0.71 




MTTB 


3.5 0 


3.51 


0.2 7 7, 


50-TH PERCNT 


2.4 3 


2.47 


1.57 % 


90-TH PERCNT 


8. 0 6 


7.25 


11.29 % 


95-TH PERCNT 


10.49 


9.83 


6.72 % 


CHI-SQR STAT 


11.20 


7.20 




DEG 0? FREED 


5 


4 




SIGNI? LEVEL 


0. 476 E-01 


3. 1 26E+00 





BOX PLOT SST41 



-I + I * 0 



30X PLOT LNSET41 



I + I 
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HISTOGRAM S2T41 



MIDDLE OF 
INTERVAL 

0 . 

2 . 

4. 

6 . 

8 . 

10 . 

12 . 

14 . 

16. 

18. 

HISTOGRAM LNSET 

MIDDLE OF 
INTERVAL 
-0.5 
0.0 
0.5 
1.0 

1.5 

2.0 

2.5 
3.0 



OF 

ATIONS 

***** 

* * * * * ***** *** 
***** ***** 
**** 

* 

* 



* 



OF 

ATIONS 

*** 

**** * ** 
***** ** 
*** 

***** **** 
***** 

* 



NUMBER 

OBSERV 

5 

13 

10 

4 

1 

1 

0 

0 

0 

1 

41 

NUMBER 

OBSERV 

3 

7 

7 

3 

9 

5 

0 

1 
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‘fcfcpAia Tir-*** 
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APPENDIX B 
T able 6 

TOTAL RESULTS (46 cases) 

Exponential Distribution : Only Chi Square test used. 

Accept Reject 

6 | 40 i 

I I 



Lognormal Distribution : Chi Square, Kolmogorov Smirnov 

and W test used. 



Chi Sg. 
test 







K-S 


test 








A 




R 




w t 


es t 


w 


test 




A 


R 


A 


R 


A 


1 12 


1 


7 


i ! 




1 






1 


R 


! i 


2 


3 


7 1 




1 






1 



( 34 cas es ) 



No W test 



K-S test 
A R 



Ch. So. 
test 



A| 2 | 2 | 

a! 0 j 8 | 



( 12 cases ) 



Gamma Distribution : Kolmogorov Smirnov test used. 



H.O.H. estimators 
A R 







l 


i 


M. L. S. 


A 


19 | 

1 


3 1 

1 


estimators 


a i 

i 


! o ! 

i i 


19 ! 
1 
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Table 7 

POMP DATA RESULTS (7 sets) 

Exponential Distribution : Only Chi Square test used. 

Accept Reject 

s ri~ i 

l i i 

Lognormal Distribution : Chi Square, Kolmogorov Smirnov 

and W test used. 











^ rz o v 








M r J 




u .est 






w 


A 

tes t 


w 


R 

test 






i* O « 








A 


R 


A 


R 






K-S 

n 


test 




A 


! o 
1 


1 


2 


o i 
1 






A 






Chi Sq. 


Ch. So. 
test 


A 


0 


0 


test 


R 


1 o 

1 


1 


0 


2 1 
1 


R 


0 


1 




















( 6 


cas es 


> 






( i 


case ) 



Gamma Distribution : Kolmogorov Smirnov test used. 



M.O.M. estimators 
A R 



W.L.E. A | 3 

1 


1 


I 

estimators R j 0 | 

1 1 


1 6 ! 

1 1 
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I able 8 

REPORT DATA RESULTS (13 sets) 

Exponential Distribution : Only Chi Square test used, 

Accept Reject 

I 1 S 12 i 

! I I 



Lognormal Distribution : Chi Square, Kolmogorov Smirnov 

and W test used. 



K-S test 



Chi Sq. 
test 



A 

W t es 


4. 

u 


w 


R 

test 


A 


R 


A 


R 


! - 
1 


0 


2 


o ! 

i 


i 0 


0 


2 


3 


( 


11 


cas es 


> 



No W test 



K-S test 
A R 



Ch. Sq . 
test 



1 

1 - 


0 1 


i i i 


1 

1 


o ! 


i i 

i i 



( 2 cases ) 

Gamma Distribution : Kolmoaorov Smirnov test used. 



N.O.M. estimators 
A R 

M.L.E. A | 7 

!_ 

• „ 

estimators R | j | 3 
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Tab Is 9 

FIELD DATA RESULTS (26 sets) 



Exponential Distribution : Only Chi Square test used. 

Accept Reject 

I 

I 4 
I 




Lognormal Distribution : Chi Square, Kolmogorov Smirnov 

and » test used. 



Chi Sq. 
test 





K-S 


te st 


w 


A 

test 


R 

W test 


A 


R 


A R 


i a 

i 


0 


3 1 ! 

1 


i 


1 


1 2 1 

1 



( 17 cases ) 



Ch. Sq 
test 



Gamma Distribution : Kolmogorov Smirnov 



No W test 



K-S test 


A 


R 


A | 2 


1 


| — — — — — 

R | 0 


6 



( 9 cases ) 
test used. 



M.O.a. as ti manors 
A R 



ft. L. E. 


A 


11 1 
1 


5 ! 

I 1 


estimators 


R 


o ! 
1 


! 

1 



S7 



sample sizes) 



Table 10 

TOTAL % RESULTS (w.r.t. 



For exponential distribution : 



Sample Size 


A 


3 


20-29 


30. 8 


59.2 


30-39 


3.3 


9 1.7 


40-50 


0 . 


100. 


GT. 50 


8.3 


91.7 



For lognormal distribution : 



Sample 

Size 


Chi- Sg. test 
A R 


W 

A 


test 

R 


K-S 

A 


* 4 - 

* 


20-29 


53.8 46.2 


71.4 


28.6 


46. 2 


53.9 


30-39 


66.7 33.3 


58.3 


41.7 


50. 


50. 


40-50 


58. 3 4 1.7 


58.3 


41.7 


33.3 


66.7 


GT. 50 


25. 75. 


— 


— 


16.7 


83.3 


gamma distr 


ibutior. : 










Sample 


A in 


H in 


A in 


MLE 




Size 


both 


both 


R in 


MOM 




20-29 


84.6 


7.7 


7. 


7 




30-39 


33.3 


50. 


16. 


7 




40-50 


25. 


52.5 


12. 


5 




GT. 50 


16.7 


15.7 


66. 


6 
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APPENDIX C 



COMPUTER PROGRAMS FOR DATA ANALYSIS 

1. LGHBHL (FORTRAN) [Ref. 1] 

The program makes a Chi-Squared Goo dress -Of- Fit test for 
the exponential and lognormal assumptions. The main program 
computes expected values of ordered observations, sample 
mean and standard deviation based on the maximum likelihood 
estimates of the parameters of the lognormal distribution, 
calculates the 50th, 90th, 95th percentiles of the exponen- 
tial and lognormal distribution assumptions, percentage 
differences between the calculated parameters of the 
exponential and lognormal distributions. 

2. BOX PLOT (MINI TAB) [Ref. 5] 

Syntax: Boxplot Vector 

Description: Generates a boxplot display for a vector 

of data. A rectangular box with ends corresponding to lower 
and upper quartiles is presented with the median marked with 
an plus sign, and lines are drawn on each side of the box 
with crosses marking the lowest and highest data values wit- 
hin half guartile distance in each side. Outliers are mark- 
ed with stars and circles. 

3. NORMPLOT (APL) [Ref. 21 ] 

Syntax: Normplot vector 

Parameters: Wid- controls the horizontal size of display 
Dep- controls the vertical size of display 
Subprograms: NSCORSS, PLOT 

Description: Generates a normal probability plot for a 

vector of data. If data fits normal then plot is straight 
lin e. 

4. EXPONPLOT (APL) [Ref. 21] 

Syntax: Sxponplot vector 

Parameters: Wid- controls the horizontal size of display 
Dep- controls the vertical size of display 
Subprogram: PLOT 
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Description: Generates a exponential probability plot 

for a vector of data. If data fits exponential then plot 
is straight line. 

5. GAM (A ?L) 

Syntax: GAM vector 

Subprograms: EST GAM , GAMMA 

Description: Perform K olmogoro v-S mirnov Good ness-of- Fit 

test and gives test statistic and '’reject' 1 / "do not reject" 
idea for gamma assumption for both MOM and MLE estimators. 

6. KOL (APL) 

Syntax: Kol vector 

Subprogram: MOM, I1JVV 

Description: Performs a Kolmogorov Smirnov Goodnees-Of- 

Fit test for a data vector for lognormal assumption and 
print test statistic. If test statistic is bigger than the 
tabulated value "reject." the assumed distribution, otherwise 
"do not reject". 

7. W (APL) 

Syntax: A W Vector 

Description: Performs a W test for the lognormal as - 

sumption. A must be a data vector (from Table IX, Ref. 2) 
and X vector that is testing against lognormality . A and X 
must have the same dimension, output will be W test statis- 
tic. Then using the value of A ,yu and 6 (Table XI, Ref. 2) 
for the appropriate sample size and than using standardized 
normal distribution determine the probability of obtaining 
value less than or egual Z, which is the significance level 
of the test. 
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<? gam can? 

7 GAM T* 

Cl 2 fl PROGRAM PERFORMS KOLMOGOROV SMIRNOF TEST FOR 

c 23 R GAMMA assumption for mom and mle estimators 
C33 ri.YC*ri 

C4] ESTGAM Y 

053 XlfYfE[l J 23 

C63 X2*-v + EC2i2J 

C7 3 n*-fY 

08 3 Z«-0*886*Hx0*5 

093 «1«-EC1J1D 

C103 «2<-EC2?n 

Cl 13 «<-0 f K*K l <N -1 ) )tN 

Cl 2 3 ifiK-0 

C 1 3 3 J<-JJ«-Q/Q 

C 1 4 3 L2;i4-x + l 

Cl 53 -»L.1X\(X>H) 

C 1 6 3 fi l gamma ;< i c i 3 
Cl 73 

C 1 3 3 -»L2 

C193 li »r«- | ( j-k ) 

C203 RR*.RCfR 3 

C213 ' test statistic for mom estimators is 1 , r P: ft c 1 3 

C223 -»L3x KRRC13 < 2 ) 

[233 ' reject for alpha=o, 05 F o p : mom estimators • 

C243 -M -5 

C253 t-3 ; * BO HOT REJECT FOR AL.RHAX0.05 FOR MOM ESTIMATORS 1 

C263 L 5 jiifii +1 
C273 ( tx >M) 

C233 a 2 gamma 2 C 1 1 3 
C2°3 JJfJJ, f 
C 30 3 -H -5 

C313 ( jj-k) 

C323 ss<-sc 9 S 3 

C333 1 test statistic for mle estimator is 1 ,tSS[i] 

C 34 3 ->L7x \ (SSC13 <2> 

C353 1 reject for alrhax0.05 f gr mle estimators • 

C 36 3 ->99 

C37] L7 : • do hot reject for alphaxQ.05 FOP: mle estimators ' 
C 33 3 499 

9 



GAM SET1 

TEST STATISTIC FOR MOM ESTIMATORS IS 0.363 
REJECT FOR ALRHA = 0.05 FOP: MOM ESTIMATORS 
TEST STATISTIC FOR MLE ESTIMATOR IS 0.154 

REJECT FOR. ALPHAxO.05 FOP ' MLE ESTIMATORS 
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C13 
C23 
C 3 3 
C 4 3 
C53 
E63 
C 7 3 
C83 

f 103 
C 1 1 3 
C 1 2 3 
Cl 33 
C 1 4 3 
C 153 
C 1 6 3 
C 1 7 3 
C 183 
C 1 9 3 
C203 
C 2 1 3 
C223 
C 23 3 
C 24 3 
C 25 3 
C26 3 
C 27 3 
C 28 3 
C 29 3 
C303 
C31 3 
C32 3 
C333 



v kol 
? kol x 



CQ3'7 



n PROGRAM PERFORMS K O L M O G O R O V - S M I ft N O V 
n 500DHE55-0F-FIT TEST FOR LOGNORMAL 
A ASSUMPTION , X MUST BE DATA ARRAY 
ORDER DATA AND FIND SAMPLE MEAN AND 

MO )■( 

A PUT DATA IN STANDARD NORMAL CASE 



2 4- ( < ) —KB ) -f S 



VAR I A N C E 



OF 



LOG DATA 



H f- f N 

F'PPfO/O 
If ilfO 

L3;iifii+i 
4 L 4 X \ ( ii > W) 

^ FIND P VALUE FOR EVERY DATA POINTS AND PUT IN ARRAY 
«-l ; 14*1 + 1 

^L4X \ ( i>?9?) 

INVV P +1x0*001 
^L2x\ (22>2[II]) 

+ *-l 

L 2 : ppp+ppp f p 

-M-3 

L4: «+< \ (N-l ) )+* 

KfOr K 

CCf ( f K ) - ( f PPP ) 

CD+(l,CC)f 1 
DCf (,CD) 

PPP+PPP , DC 

* FIND MAXIMUM ABSOLUTE DEVIATION AND PRINT AS TEST STATISTIC 



K K f- | ( P P p - K ) 

LfKK[fKK] 

'TEST STATISTIC r 1 *f L Cl] 

T v 4- Q ♦ 3 8 6 ^ ( w * 0 ♦ 5 ) 

'TABULATED VALUE - 1 , <> T 7 
V 



KOL SET 1 

TEST STATISTIC =0*146 
TABULATED VALUE =0.135 
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v w ZQ2v 

V Y W J ^ 1 ? A2 7 A 3 f s 2 7 K ? N ? NN 



CIU 


ft 


PRODUCES W TEST STATISTIC 


FOR W 


TEST 


C2H 


ft 


TrVECTOft FOR W TEST LENGTH 


OF A R G U M E N T 


C3: 


ft 


X = VECTOR THAT IS TESTING 






[4: 


ft 


X AND Y MUST BE SAME DIMENSIONS 




cs: 




N <- j> X 






C6: 




I 






Z71 


ft 


ORDER DATA 






ZQl 




X f>{[^ x 2 






Z91 


ft 


CALCULATE S*2 TERM 






cio: 




a 1 *• + / m '< 






cm 




<=>2<- + / ( •« ) *2 






C 12 : 




A 3*-< A 1 *2) 






C13H 




S2RA2-A3 






C14H 


ft 


CALCULATE B TERM 






cis: 




K 4t- l N ~ 2 






C 161 




LXf f X 






Z171 




LLXfLX[ ^LX 2 






cia: 




B «- + / Y X ( (KfLLX)-(KALX) ) 






Cl?3 


ft 


FIND W TEST STATISTIC AMD 


P R I W T 


AS Uf 


C2o: 










C2i: 




* W TEST STATISTIC = ', t Z 










A43 W SET ]_ 

w TEST STATISTIC = 0*967 
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7 GAMMA CQDV 
7 A GAMMA X j V } T J 2 | ft } O i H 

[ 1 3 ft INC, GAMMA FHC , -SHAPE=A > 0 > VECTOR AF:G , '< '■ Q , SCALE PARAM, 

C2H ft GIVES THE LEFT TAIL AREA FOR GIVEN A AMD~X 
C3D ft <=> IS SHAPE PARAMETER, X IS DATA VALUE 

C43 A*-LA 
C 53 P<-l-*-X 

C63 +oxu=a 
C 7 3 Pf(fXf»x)fO 

can vf(x<A) V x<i 
C9 3 4l-X\0 = + /V 

CIO 3 i>v/x 

[in N+r~7+*<r/T)+A+i 

C 1 2 3 Z«-l+ + /X\To 

C 1 3 3 -»< lt3 + D«-c> x x i <a 

C 1 4 3 PCV/\f X3^2 X (TRA)^( JA) X*T 
C 153 -M- 

C 1 63 F[V/lf!-!]«-(TTA)X 2 X( X / ( V X * - T -r A- 1 ) * , * A- * A- j. ) X ( Y * A- A ) * J A- A 
C 1 7 3 Lj-»0x \0 = + /"V 
C 1 8 3 r «- ( -v v ) / x 
C 1 9 3 5«-20r2x» 

C203 Z + A CF T 

C21 3 + < It3 + 0LC) X \ 1 < A 

C 22 3 F-C(" v )/\f'<3<-l-2X(f*A)X (*-T ) + : A — 1 

C 23 3 + 0 

C243 F '[('' v )/U > :<]<-l-'l'x2X(X/(Tx*-TTA-l)o,-rA-lA-l)X(TjA.A)r!A-A 

V 



V IHVV C03v 
7 IHVV A 

C13 ft 'IHVV A' RETURNS THE VALUE OF STANDARD 

C23 ft normal distribution as left tail area a 
C 3 3 & + 1 ( lxA>0,5) -A 

C43 X4-(-*&xB) *0.5 

C53 x i <- 2 . 5 1 55 17 + xx 0.802853 + 0. 01 0328 xx 
C 6 3 x 2 + 1 + x x 1 ,432783 + xxO, 1 89269 + 0 . 00 1308 x x 

C 7 3 zz«- < x i fX2 ) -x 

C83 zz<-zz-2xzzx a > o , 5 

v 



7 MOM C 0 3 7 
7 MOM X 

Cl 3 ft RETURNS sample mean and variance 

C23 XB*. ( +/'.<■) TfX 

C3 3 s<-(+/(X-XB) * 2 ) + ~ 1 + f x 

V 



104 



7 E3TGAM £ Q 3 V 
7 E5TGAM XJAjBjFjFP 

[13 fl RETURNS METH, OF MOMENTS EST, THEN MAX LIK EST, OF GAMMA 
121 fl OUTPUT AS 2x2 MATRIX its first row mom and second row mle 
131 fl estimators shape and scale parameters respectively 

[43 MOM X 

[53 A + X B B +* S + X B 

[63 Ef A , B 

171 (f/ix>Tf x 

LSI LJAQfA 

[93 F f & j- ( a A *r X B ) — o F S I & 

[10 3 F F >(tA)-1 RSI A 
£113 a+a-f+fp 
£123 ->LX\ (0.0001) < I A-*Q 

£13 3 E<- 2 2 p£,a,xb+a 

V 



v FSI £ £] 3 V 

7 F>H PS I YjCjlVJJIVJKjKKj YY J V j Z J T J I 
[13 fl N IS THE ORDER OF THE DERIVATIVE OF THE PSI FUNCTION 

[ 9 3 fl T ( > 0 ) 15 THE ARGUMENT , SCALAR OR VECTOR 

£33 c*. 10 

£ 4 3 I V 4- \ f Y f. , T 

£53 F>2fK(-( r r)po 

£63 K K + r C - Y [ J I V f ( V f f ( C ) / I v 3 

£73 4Lix \ (f IV)=P (#rV)/iv 

£83 t f t [ j i v ] 

£93 i«-0 

£103 ^jifi + i 

£113 yy<*k k [ i ]f T[ i 3 

£12 3 2[I]f(-H)XT/(( YY-! ) + l KK[I3 ) *-l+N 

£133 -M- 2 x \ x <f jiv 

£143 ZfV\z[^Jiv] 

£153 KfV\KK 
£163 l. 1 J-tSlx \N>0 
£173 Pf-(iK + Y)-( 2 xK + r) ^ — l 
£183 ->2 + 0 «-c 

£19 3 si J P4. < < • N — i ) X ( Y + K ) *-N ) + ( l N ) X0 • 5X ( Y + K ) *-n + 1 
£203 F*P( (“1 ) *N+i ) XF + 2 + N JEX Y + K 

V 



DATA 
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